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In SU(2) lattice pure gauge theory we study numerically the dependence of the ghost propagator 
G(p) on the choice of Gribov copies in Lorentz (or Landau) gauge. We find that the effect of Gribov 
copies is essential in the scaling window region, however, it tends to decrease with increasing /3. 
On the other hand, we find that at larger /3-values very strong fluctuations appear which can make 
problematic the calculation of the ghost propagator. 



I. INTRODUCTION 

The nonperturbative study of the ghost propagator is 
of great interest for the understanding of the mechanism 
of confinement. The Kugo-Ojima confinement critcrium 
Q, Q] is formulated in terms of the ghost propagator 
G ab {p) at p — > and expresses the absence of colored 
states from the spectrum of physical states. Moreover, 
confinement of gluons has been related more directly 
to the suppression of the gluon propagator in the limit 
p — > 3]. In both cases, the propagator in question is 
defined in the Landau (or Lorentz) gauge. 

In a series of papers Zwanziger |^ has suggested that 
this behavior might result from the restriction of the 
fields in the transversal plane r = {^4 : d^A^ = 0} to 
the Gribov region = {A : d^A^ — 0, M > 0}, where 
the Faddeev-Popov operator M is non-negative. 

From studies of the coupled Dyson-Schwinger equa- 
tions for gluons and ghosts 0, Q] it is well-known that 
the infrared behaviour of gluon and ghost propagators 
is closely related pj: the gluon propagator D^ b (p) = 
{S^v —PuPv/p 2 ) Z g i(p 2 )/p 2 is damped in the infrared like 
Z g i(p 2 ) oc (p 2 ) 2K , while the ghost propagator G ab (p) — 
S ab G(p) — 5 ab Z g h(p 2 )/p 2 is more singular than the free 
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propagator, Z g h{p 2 ) oc (p 2 ) K . In a particular truncation 
scheme K — 0.595 has been determined 

There are only relatively few previous lattice stud- 
ies of the ghost propagator Q, 0, in con- 
trast to numerous investigations of the gluon propaga- 
tor As for the latter, is not yet clear from 
the lattice Jjj whether Z g i(p 2 )/p 2 -> or ^ with 
p 2 — ► The lattice volumes might still be insuf- 
ficient to decide this question. The singular behavior of 
Z„h(p 2 ) is seen to become stronger with increasing vol- 
ume pjj. This supports the expectation 0] that the 
sample of physically important gauge field configurations 
A G r, which constitutes the Euclidean functional inte- 
gral, in the thermodynamical limit V — > oo is concen- 
trating towards the edge of the Gribov region, the first 
Gribov horizon <9il where the lowest non- vanishing eigen- 
value of the Faddeev-Popov operator is approaching zero. 
This statement is the content of Zwanziger's horizon con- 
dition 0, which can be related to the Kugo-Ojima 
criterion. 

All this is complicated by the non-uniqueness, first 
pointed out by Gribov 3], of the intersection with T of 
the gauge orbit A 9 of any gauge field A, even if restricted 
to the Gribov region fi. Practically, the Landau gauge 
is implemented by maximizing (with respect to gauge 
transformations g) a certain gauge functional. Usually, 
such a problem leads to more than a single maximum, 
which are gauge copies (Gribov copies) of each other, 
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hence to a non-unique definition of gauge dependent ob- 
servables. Thus, in a lattice investigation one has to de- 
termine which observables are really subject to the so- 
called Gribov problem which reflects the dependence of 
an observable on the restriction (if possible) to the copy 
corresponding to the absolute maximum of the gauge 
functional. More precisely, one has to study whether 
this dependence disappears when one is approaching the 
continuum and/or infinite volume limit. Otherwise this 
would indicate the persistence of a real Gribov problem 
to which Gribov has drawn the attention. On the lattice, 
the structure of the Gribov region has been closer investi- 
gated under this aspect only by Cucchieri 21] some years 
ago. 

Here we are mainly dealing with the infrared behav- 
ior of the calculated ghost propagator. In the result of 
a study for SU(2) gluodyamics |llj, Cucchieri came to 
the conclusion that the ghost propagator depends on the 
selection of the highest among more and more maxima 
of the gauge functional while the gluon propagator does 
not depend. This study was restricted on one hand to 
the strong coupling region (/3 = 0.0, 0.8, 1.6) where these 
observations apply, and /3 = 2.7 where no gauge copy de- 
pendence was seen at all. These (3 values are outside the 
physically interesting scaling region. In a more recent 
paper 13], it has been reported that the gauge copy de- 
pendence of the ghost propagator in the more interesting 
scaling region (at (3 = 2.15, 2.2, 2.3 and 2.4 for lattices 
16 3 x 32) has been found to be within the statistical er- 
rors, on a level which is called Gribov noise. 

In the present paper we reanalyse the scaling region 
at f3 = 2.2, 2.3, 2.4, 2.5 and 2.6 for lattices 8 4 and 16 4 
by comparing two ensembles of gauge-fixed field config- 
urations. One ensemble ("fc") consists of an arbitrary 
maximum (usually the first being found), and the other 
consists of the best (relative) maximum ("be") among 
N copy local maxima of the gauge functional. We find 
that the difference of the ensemble averages of the ghost 
propagator for the lowest non-vanishing lattice momen- 
tum between the two ensembles does not vanish, except 
for the highest (3 value. Hence the Gribov problem re- 
mains a serious obstacle for a unique definition of the 
SU(2) ghost propagator in the scaling region. More seri- 



ous is an unexpected observation in the higher-/? region. 
We find an intermittent behavior of the ghost propagator 
estimator for the lowest non-vanishing momentum, sig- 
nalled by anomalously large, isolated fluctuations of the 
ghost propagator G(p m i n ) (see below) within the time 
history of uncorrelated configurations. We stress already 
here that this behavior is not a Gribov copy problem 
since the anomalous peaks of G{p m i n ) are observed both 
for the first and the best Gribov copy, entering the ensem- 
bles "fc" and "be", respectively. We have tested whether 
this is correlated with various infrared observables. For 
the time being, two hypothetic causes must be excluded 
as a viable explanation of the phenomenon. 

In Section 2 we recall the definition of the gluon field 
A^, the definition of the Lorentz (or Landau) gauge, the 
structure of the Faddeev-Popov operator and the defini- 
tion of the ghost propagator. Details of the simulations, 
the gauge fixing and the observation of Gribov copies are 
reported in Section 3. In Section 4 we discuss the results 
on the ghost propagator. We conclude in Section 5. 



II. FADDEEV-POPOV OPERATOR AND 
GHOST PROPAGATOR 

A. Definition of the gluon field and Faddeev-Popov 
operator 



For the Monte Carlo generation of ensembles of non- 
gauge-fixed gauge field configurations we use the stan- 
dard Wilson action [22|, which for the case of an SU(N) 
gauge group is written 



S 



I 



X fl~>V 



1 

N 



Re Tr 



(^u x ^u x -^ fJ _- L ,ul +l/ . tl u^ 1 ^j 



(1) 



Here go is a bare coupling constant and U xfl G SU(N) 
are the link variables. The field variables U Xfl transform 
as follows under gauge transformations g x : 



U 



Xfl. 



U L = g x U xti g x+h 



g x eSU(N). (2) 



For SU(2) gauge links U xfJ _, a standard definition 23] of 
the lattice gauge field (vector potential) ^4 x +/i/2, A i is 



A 



x-\-p,/2,fi 



(3) 
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Therefore, for SU(2), the link can be written 

U xl i = b xfl 1 + i b xf , a = b xfl 1 + i At+aAm 



(4) 



In lattice gauge theory the usual choice of the Landau 
gauge condition is Q 

4 

[8A) X = ^ (Ar+AAp - Ax-fi/2;v) = , (5) 

which is equivalent to finding an extremum of the gauge 
functional 



(6) 



with respect to gauge transformations . After replac- 
ing U => U 9 at the extremum the gauge condition (J5J 
is satisfied. In what follows this gauge is referred to as 
Landau gauge. 

The lattice expression of the Faddeev-Popov operator 
M ab corresponding to M ab — —d^D°f in the continuum 
theory (where D ab is the covariant derivative in the ad- 
joint representation) is given by 

M± = 



ab _ f / qab , qab \ r 

\ u x)j. xfj.) u y\x+ii 
~ {^x-fi;fj. + -^-x-fi-.fj) $%r,X—(i | 



(7) 



where 



M 9 ^:r+/i/2;M 



cafc _ rafc _ m r f . xai, 



(8) 



From the form JSJ it follows that a trivial zero eigen- 
value is always present, such that at the Gribov horizon 
dT the first non-trivial zero eigenvalue appears. Con- 
versely, it is easy to see that for constant field configura- 
tions, with b xfl = 6° and b xfl = b® independent of x, there 
exist eigenmodes of M with a vanishing eigenvalue. Thus, 
if the Landau gauge is properly implemented, M[C7] is a 
symmetric and semi-positive definite matrix. 

B. Ghost propagator 

The ghost propagator G ab (x,y) is defined as 0,0 

G ab (x, y) = 8 ab G(x-y) = {( )»* [U] ) , (9) 



where M[U] is the Faddeev-Popov operator. Note that 
the ghost propagator becomes translational invariant 
(i.e., dependent only on x — y) and diagonal in color 
space only in the result of averaging over the ensem- 
ble of gauge-fixed representants (first or best gauge-fixed 
copies) of the original Monte Carlo gauge configurations. 
The ghost propagator in momentum space can be written 
as 

G{ V) = ^J2 e ~ 2niP<X ~ V) (( M ~ X )l a yM) > ( 10 ) 
V 

where V — L 4 is the lattice volume, and the coefficient 
7^7 is taken for a full normalization, including the indi- 
cated color average over a = 1, .., 3. 

We mentioned above that M[Z7] is a symmetric and 
semi-positive definite matrix. In particular, this matrix 
is positive-definite in the subspace orthogonal to constant 
vectors. The latter are zero modes of M[U}. Therefore, 
it can be inverted by using a conjugate-gradient method, 
provided that both the source ip a (y) and the initial guess 
of the solution are orthogonal to zero modes. As the 
source we adopted the one proposed by Cucchieri 11]: 

r(y) = s ac e 2 ^-y p=/= (0,0, 0,0), (11) 

for which the condition ~^2 y ip a (y) = is automatically 
imposed. Choosing the source in this way allows to save 
computer time since, instead of the summation over x 
and y in Eq. I|10[l . only the scalar product of M -1 ^ with 
the source tp itself has to be evaluated. In general, the 
gauge fixed configurations can be used in a more efficient 
way when the inversion of M is done on sources for c — 
1,..,3 such that the (adjoint) color averaging, formally 
required in Eq. I|1(J|) , will be explicitely performed. 

III. SIMULATION DETAILS 

The numerical simulations have been done for SU (2) 
pure gauge theory using the standard Wilson action, for 
lattice volumes L 4 with L = 8 and L — 16. At a given 
lattice size L for each j3 value we have generated N con f 
independent mother configurations, for which the Lan- 
dau gauge was fixed N copy = 20 times, each time start- 
ing from a random gauge transformation of the mother 
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configuration, obtaining in this way N copy Landau-gauge 
fixed copies. 

Two consecutive configurations (considered as inde- 
pendent) were separated by 100 and 200 sweeps for lat- 
tice sizes 8 4 and 16 4 , respectively. Each sweep consisted 
of one local heatbath update followed by 4 or 8 micro- 
canonical updates for 8 4 or 16 4 lattices. In all our 
runs we have measured the integrated autocorrelation 
time for the plaquette, for the Polyakov loop and for the 
ghost propagator (separately for each momentum p) . In 
all cases, the relation r int ~ 0.5 was observed, showing 
that the consecutive configurations are effectively inde- 
pendent. 

The actual measurements of the ghost propagator were 
done for the "first", i.e. in fact an arbitrary gauge copy 
and for the "best" one among the N copy copies. If the 
first copy turned out to be the best, the ghost propagator 
was measured only once, and the result simultaneously 
entered the two different gauge-fixed ensemble averages. 
In the following the two ensembles are labelled "fc" and 
"be", refering to the first or the best gauge copy, re- 
spectively. In Tabled we give, for each set of simulation 
parameters (L,/3), the number of times the first copy 
produced turned not out to be the best, i.e., did not cor- 
respond to the relative maximum of Fjj(gi) among the 
N copy copies. 

As the gauge fixing procedure we used standard Los 
Alamos type overrelaxation with u> = 1.7. The iterations 
have been stopped when the following transversality con- 
dition was satisfied: 



4 

max 

x, a I — ' 



A a 



A a 



(12) 



We used the parameters ei or — 10~ 10 or 10~ 9 for lat- 
tice size 8 4 or 16 4 , respectively. In our test runs it was 
found that further decreasing e; or does not affect the re- 
sults for the ghost propagator. Also it was checked that 
these values of ei or are sufficient for identifying, accord- 
ing only to the values of Fjj(gi), Gribov copies which are 
actually global gauge transformations of each other and 
conversely for distinguishing this from the case of actu- 
ally inequivalent lattice Gribov copies. 

In Tabled for each set of simulation parameters (L, /3), 



size 


(3 


prconf IM con ! 
multiple copies/ 


* j-copies I jtj-copies 
nonequiv/ total 




8 4 


1.6 


500 / 500 


8263 / 10000 


446 




2.0 


490 / 500 


4431 / 10000 


354 




2.1 


468 / 500 


3460 / 10000 


311 




2.2 


426 / 500 


2180 / 10000 


235 




2.3 


301 / 500 


1364 / 10000 


150 




2.4 


184 / 500 


877 / 10000 


92 


16 4 


2.0 
2.1 
2.2 
2.3 
2.4 
2.5 


25 / 25 
25 / 25 
25 / 25 
25 / 25 
25 / 25 
24 / 25 


500 / 500 
500 / 500 
500 / 500 
494 / 500 
337 / 500 
169 / 500 


25 
25 
23 
25 
23 
14 



TABLE I: The Table shows in the 3rd column the num- 
ber N mu/ti P ie copies of configurations, for which non- equivalent 
Gribov copies have actually been obtained, out of a total num- 
ber N con f of configurations which underwent gauge fixing; in 
the 4th column the total number N^°J^^ uiv of non-equivalent 
Gribov copies out of a total number N^°^ B = N con f x N CO py 
of gauge copies under investigation. The last column presents 
the number of times out of N con ^ that the first copy was not 
identical to the best (relative maximum) copy. 



we present also the number of configurations for which 
Gribov copies have been found and the total number of 
different Gribov copies. 

The momenta p for the propagator G(p) were taken 
with all spatial components put equal to zero: p = 
(0, 0, 0, ki/L), where £4 was restricted to k± — 1, 2, 3, 4. 

IV. DISCUSSION OF THE RESULTS 

From Table UJ one can learn that at the lattice size 
8 4 the fraction of Monte Carlo configurations which are 
represented by more than one gauge-fixed configurations 
(among 20 attempts to find copies) drastically begins to 
decrease at (3 = 2.3. Parallel to this also the multiplicity 
of actually different copies among 20 drops down. The 
decrease of the number of available basins of attraction 
for the gauge fixing process is a finite- volume effect. For 
the bigger lattice size (16 4 ) one sees that the fraction of 
Monte Carlo configurations with more than one gauge- 
fixed configurations practically does not depend on (3. 
However, the multiplicity of non-equivalent copies among 
the 20 obtained copies starts to decrease from (3 = 2.3. 

From Table ITU one can see for separate small momenta, 



8 4 lattice 



p 




Copy 


fc 4 = 1 


k 4 = 2 


&4 = 3 


k 4 = 4 


1.6 


500 


be 


6.58(4) 


1.327(5) 


0.628(2) 


0.501(1) 


1.6 


500 


fc 


7.02(6) 


1.363(5) 


0.638(1) 


0.508(1) 


2.0 


500 


be 


5.15(3) 


1.013(2) 


0.491(1) 


0.3970(4) 


2.0 


500 


fc 


5.46(9) 


1.028(3) 


0.495(1) 


0.3995(6) 


2.1 


500 


be 


4.62(3) 


0.920(2) 


0.4545(6) 


0.3701(4) 


2.1 


500 


fc 


4.89(7) 


0.935(3) 


0.4573(8) 


0.3719(5) 


2.2 


500 


be 


4.06(3) 


0.823(2) 


0.4189(4) 


0.3444(3) 


2.2 


500 


fc 


4.26(4) 


0.833(2) 


0.4203(5) 


0.3450(3) 


2.3 


500 


be 


3.60(4) 


0.744(1) 


0.3903(3) 


0.3238(2) 


2.3 


500 


fc 


3.65(4) 


0.747(2) 


0.3909(4) 


0.3241(2) 


2.4 


500 


be 


3.38(5) 


0.691(1) 


0.3710(4) 


0.3098(2) 


2.4 


500 


fc 


3.47(7) 


0.692(2) 


0.3712(4) 


0.3099(2) 


16 4 lattice 


P 


i<>meas 


Copy 


fc 4 = 1 


k 4 = 2 


&4 = 3 


k 4 = 4 


2.2 


296 


be 


20.1(1) 


3.87(1) 


1.494(2) 


0.8078(6) 


2.2 


296 


fc 


21.3(1) 


3.97(1) 


1.509(2) 


0.8115(6) 


2.3 


270 


be 


17.3(1) 


3.29(1) 


1.303(1) 


0.7248(4) 


2.3 


270 


fc 


18.0(1) 


3.33(1) 


1.310(2) 


0.7268(5) 


2.4 


370 


be 


14.8(1) 


2.83(1) 


1.165(1) 


0.6673(3) 


2.4 


370 


fc 


15.6(1) 


2.87(1) 


1.171(1) 


0.6690(3) 


2.5 


294 


be 


13.7(2) 


2.56(1) 


1.088(1) 


0.6353(3) 


2.5 


294 


fc 


13.9(2) 


2.58(1) 


1.090(1) 


0.6358(3) 


2.6 


229 


be 


13.6(4) 


2.41(1) 


1.043(2) 


0.6161(5) 


2.6 


229 


fc 


13.8(4) 


2.41(1) 


1.044(2) 


0.6164(5) 



TABLE II: The ghost propagator G(p) from Eq. as 
a function of k 4 — 1,2,3,4. We have set the momentum 
p — (0, 0, 0, k 4 /L), where L = 8, 16 is the lattice size. The 
averages over the gauge configurations in Eq. (I10H were taken 
in two different ways: "fc" means the average taking only the 
gauge-fixed copy generated first for each configuration, "be" 
means the average over only the best (relative maximum) copy 
among 20 different gauge-fixed copies that we have generated. 

how the average value of the ghost propagator differs be- 
tween the two ways to deal with the Gribov copy prob- 
lem: to ignore it (N copy — 1) or to inspect N copy = 20 
copies. For all momenta, the ensemble consisting of the 
first copies ('fc') turns out to give slightly larger values 
than the ensemble including always the best copy ('be'). 
For the lowest non- vanishing momentum this is shown in 
Fig. [l]for the lattice 8 4 . It is visible that for j3 G [1.5, 2.4] 
the difference of G(p m j„) between the two ways of aver- 
aging is clearly outside the statistical error. 

In Fig. |21 for the bigger lattice 16 4 , the ghost propa- 
gator values for the two lowest momenta are compared 



with respect to the dependence on Gribov copies for 
(3 £ [2.2, 2.6]. Whereas for the lowest momentum the re- 
sults resemble those of the smaller lattice, for the second 
lowest momentum they are practically indistinguishable 
at the given scale. For increasing j3 the difference be- 
comes of the order of the statistical error (Gribov noise). 
At (3 = 2.6 the ghost propagator data even for the lowest 
momentum fall together within error bars. This indicates 
that the Gribov problem has disappeared for the ghost 
propagator there. 

Instead, at (3 = 2.6 a new problem arises which can 
be recognized already in Fig. [2] where we also demon- 
strate how, at (3 = 2.6, the average for the ghost prop- 
agator at the lowest momentum would be influenced by 
the removal of " exceptional configurations" . These are 
signalled as spikes in the Monte Carlo time histories of 
the corresponding observable shown in Fig. |2]for (3 = 2.6. 
Precursors of this phenomenon are visible there at lower 
(3, too, but for (3 — 2.6 the effect becomes notable. We 
notice that these spikes occur in the first as well as in the 
best gauge-fixed copy. Therefore, the existence of these 
"exceptional configurations" is definitely not a result of 
gauge fixing. 

In order to explore what the essence of these "excep- 
tional configurations" is, we have looked for correlations 
with certain "toron" excitations on one hand and with 
different Polyakov loops on the other. 

In the first case we followed the procedure applied by 
Kovacs [25I for extracting the toron content of Monte 
Carlo gauge field configurations |2j|. We evaluated for 
all four directions fj, on the lattice the corresponding 
holonomies over a /x-slice fixed at = 1 

L-l 

M x ) = n u *+*w ■ (!) 

s=Q 

We averaged this quantity over the /z-slice, 

^ = zs E (2) 

These gauge dependent quantities were normalized to 
SU(2) in the usual way 

V p P^/^detP^ . (3) 
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FIG. 1: The (i dependence and Gribov copy sensitivity of the ghost propagator G(p m in) at minimal momentum p„ 
8 4 lattice. Filled symbols correspond to the 'be' ensemble, open symbols to the 'fc' ensemble (see the text). 
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FIG. 2: The f3 dependence and Gribov copy sensitivity of the ghost propagator G(p) at minimal momentum p m in and next- 
to-minimum momentum 2p m i„ on the 16 4 lattice. Filled symbols correspond to the 'be' ensemble, open symbols to the 'fc' 
ensemble. At (3 = 2.6 the minimal momentum ghost propagator is sensitive to the inclusion or exclusion of "exceptional 
configurations" (see text). 



Then the anticipated homogeneous toron field is given 
by links U xtl independent of x, which are required to 
reproduce as follows : 

(U X ,) L = V, . (4) 



The corresponding toron gluon field can be extracted as 
We have plotted the time history of the lowest momen- 
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FIG. 3: Monte Carlo time histories of the ghost propagator G(p m in) for various /3 on the 16 4 lattice. The frequency of the 
occurrence of "exceptional configurations" increases with higher f3. 



turn ghost propagator together with the toron observable 



of the average Polyakov lines L„ for fx = 4 (middle) and 
fi = 1 (lower panel). 



T — Tr ( A toron , ) 2 



(6) 



o=l 



defined separately for the four Euclidean directions. We 
noticed that the previously mentioned spikes ("excep- 
tional configurations" ) occur independent of spikes of this 
toron observable in each of the Euclidean directions. We 
demonstrate this in Fig. 0] which shows the Monte Carlo 
history of the lowest-momentum ghost propagator (up- 
per panel) together with the histories of the toron fields 
for [i = 4 (middle) and fi = 1 (lower panel). 
We also checked the Monte Carlo sample for eventual 
correlations with the average Polyakov loop 



W = 2T1- ^ 



(7) 



Similarly, we illustrate in Fig. [5] that there are no cor- 
relations between the spikes of the lowest momentum 
ghost propagator with extremal fluctuations of the av- 
erage Polyakov loop in any of the four directions. Shown 
in the Fig. are, beside the history of the lowest- 
momentum ghost propagator (upper panel), the histories 



V. CONCLUSIONS 

In this work we studied numerically the dependence of 
the ghost propagator G{p) in pure gauge 5/7(2) theory 
on the choice of Gribov copies in Lorentz (or Landau) 
gauge with the special focus on the physically interesting 
scaling region. All simulations have been performed on 
the 8 4 and 16 4 lattices. 

We found that the effect of Gribov copies is essential in 
the scaling window region. Therefore, the Gribov prob- 
lem remains a serious obstacle for a unique definition of 
the SU(2) ghost propagator in the scaling region. How- 
ever, it tends to decrease with increasing f3 values. 

Another - and more serious - problem is presented 
by the unexpected observation, in the higher-/? region, 
of anomalously large, isolated fluctuations of the ghost 
propagator G(p m i n ) within the time history of uncor- 
related configurations. These strong fluctuations make 
problematic the calculation of the ghost propagator. 

We believe that this problem deserves a more thorough 
study, in particular how to interpret the relevant config- 
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FIG. 4: Monte Carlo time histories of the ghost propagator G(p m in) (above) for (i — 2.6 on the 16 4 lattice, compared with 
the histories of the toron fields T4 (middle) and T\ (below). The fluctuations of the latter have been arbitrarily rescaled and 
shifted for better visual inspection. 

150 i , 1 , 1 , 1 , 1 , 1 , 1 




-100 1 1 1 1 1 1 1 1 1 1 1 1 1 
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FIG. 5: Monte Carlo time histories of the ghost propagator G(p m in) (above) for [3 = 2.6 on the 16 4 lattice, compared with the 
histories of the average Polyakov lines L4 (middle) and L\ (below). The fluctuations of the latter have been arbitrarily rescaled 
and shifted for better visual inspection. 
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